A Spring-block Model for Barkhausen Noise 
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A simple mechanical spring-block model is introduced for studying magnetization phenomena and 
in particularly the Barkhausen noise. The model captures and reproduces the accepted microscopic 
picture of domain wall movement and pinning. Computer simulations suggest that this model 
is able to reproduce the main characteristics of hysteresis loops and Barkhausen jumps. In the 
thermodynamic limit the statistics of the obtained Barkhausen jumps follows several scaling laws, 
in qualitative agreement with the experimental results. The simplicity of the model and the invoked 
mechanical analogies makes it attractive for computer simulations and pedagogical purposes. 
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I. INTRODUCTION 

Barkhausen noise (BN) belongs to the family of the 
so-called crackling noises [1]. It appears as a series of 
discrete and abrupt jumps in the magnetization when 
a ferromagnetic sample is placed under varying external 
magnetic field. It is believed that the BN is a consequence 
of the fast movement of domain walls between pinning 
centers, which are either defects or impurities in the ferro- 
magnetic sample. The present paper intends to introduce 
a simple and successful mechanical model for describing 
this classic magnetization phenomenon. The simplicity 
of the model and the invoked mechanical analogies make 
this model attractive for computer simulation studies and 
pedagogical purposes. 

The Barkhausen phenomenon is interesting from sev- 
eral points of view. From a practical side, by measur- 
ing the BN there is a possibility for non-destructive and 
non-invasive material testing and control. On the other 
hand, from a pure conceptual viewpoint, by studying the 
BN one might reach a better understanding of the com- 
plex dynamics of domain walls during magnetization pro- 
cesses. 

Since its discovery (1917) BN has been intensively 
studied. Numerous measurements were done to clarify 
the statistical properties of the BN [2-4] . Regarding the 
nature of the Barkhausen noise (white noise, 1// noise, or 
l// 2 noise) even the experimental results are incoherent 
with each other. The most extensive measurements and 
data analysis were performed by Spasojevic et al. [2] with 
a commercial VITRO VAC 6025-X metal glass (quasi 2D) 
sample. After performing the statistical and numerical 
calculations, they have found: (i) power-law type distri- 
bution for signal duration with scaling exponent —2.22, 
(ii) power-law behavior for signal area with scaling ex- 
ponent — 1.77, and (iii) power-law type power spectrum 
with scaling exponent —1.6 to —1.7. From here they 
concluded that BN is not pure 1//, nor l// 2 (Brownian) 



type noise, but something between these two. Plewka, et 
al. [3] performed measurements (and calculations) with 
a similar experimental setup on an amorphous ribbon 
in an open magnetic circuit. They obtained instead a 
value around —0.9 for the scaling exponent of the power 
spectrum. From this result they concluded that BN is 
typically 1// noise. O'Brien and Weissman [4] performed 
measurements with a SQUID magnetometer on an amor- 
phous iron-based metallic alloy (2605TCA) and they sug- 
gested that BN is much closer to a white noise than a 1// 
noise and differs sharply from most typical 1// noises. 

BN received a special attention in the context of self- 
organized criticality. Self -organized criticality (SOC) is 
a term used for a class of complex phenomena where 
non-equilibrium broadband noise in driven systems re- 
flects a type of self-organization, producing states with 
power-law correlations closely analogous to critical phe- 
nomena [5]. Some of the ingredients of SOC were known 
to be potentially relevant to BN. In some cases, mag- 
netization changes have been directly observed to occur 
via avalanche process in the domain topology [6]. These 
avalanches exhibit some scaling effects, at least over a 
narrow range of parameters, and their behavior has been 
described by a SOC model [7,8]. There are however other 
approaches that put under doubt the relevance of the 
SOC concept to BN. O'Brien and Weissman [4] for ex- 
ample argued that the presumed 1/ / nature of BN and 
the observed power-law distributions are not necessar- 
ily evidences of SOC, but rather the consequences of the 
scaling properties of quenched disorder in the material. 

Many conceptually different models were elaborated 
to explain BN and its scaling properties. Without the 
intention of making a complete review, here we will men- 
tion only a few selected theoretical approaches. Alessan- 
dro et al. [9] proposed a single degree of freedom model 
(ABBM model) that considers the motion of a single do- 
main wall in a spatially rough coercive field created by 
the defects. They concluded that a mean-field approxi- 
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mation is adequate, and found power-law behavior for the 
Barkhausen pulse size distribution. Another model [13] 
which is strongly related to the previous one considers the 
motion of a single flexible domain wall in an uncorrelated 
disordered medium. This approach leads to a power-law 
distribution of the avalanche sizes (exponent —1.5) and 
durations (exponent —2), and yields an exponent —2 for 
the scaling of the power spectrum. Perkovic, Dahmen 
and Sethna [11] described the BN in terms of avalanches 
near a critical point. They used the zero-temperature 
random field Ising model (RFIM), in which the effect of 
the pinning centers was taken into account as a normally 
distributed local random field. This model was able to 
account for the power-laws characteristic for the distri- 
bution of avalanche sizes, signal area and signal duration. 
Another theoretical attempt by O. Narayan [12] consid- 
ers a multiple degrees of freedom model, studying the 
relaxational dynamics of a single domain wall in a two- 
dimensional Ising system. This model yields a power-law 
with critical exponent —1.5 for the power spectrum. The 
model predicts that in other dimensions different critical 
exponents are expected, for example in one-dimension 
the critical exponent for the power spectrum should be 
zero. This result is in contradiction with the prediction of 
the mean-field approximation for the single degree of free- 
dom model [9], which predicts the value —2 for the expo- 
nent, independently of the dimensionality of the model. 

Despite of the numerous complex models and concep- 
tually different approaches, presently none of them is able 
to face realistically all the experimentally observed fea- 
tures of the Barkhausen phenomenon. The spring-block 
model which will be introduced in this paper offers an 
elegant alternative. As it will be shown in the following, 
this model has the potential to explain within a relatively 
simple manner most of the experimentally observed scal- 
ing laws characteristic to BN. 

II. THE SPRING-BLOCK MODEL 

The model is essentially a one-dimensional spring- 
block system. It is aimed to reproduce the accepted mi- 
croscopic picture of domain wall dynamics for 180 degree 
Bloch- walls which separate inversely oriented (+| — | + 
| — | + . . .) magnetic domains (Fig.l). 

We assume that the domain walls are pinned by de- 
fects and impurities, and cannot move unless the resul- 
tant force acting on them is bigger than the strength of 
the F p pinning force. When the resulting force is greater 
than the pinning force, the wall simply jumps in the re- 
sulting force direction on the next pinning center. Apart 
of this pinning force there are two other types of forces 
acting on each domain wall. To understand these forces 
let us consider the i-th wall (which separates the (i — l)-th 
and i-th domain) free to move and all other walls fixed. 
One of the forces acting on the domain wall, Fh , results 
from the magnetic energy of the domains i and (i — 1) 



in an external magnetic field. Let us consider the exter- 
nal magnetic field as sketched in Fig.l. The interaction 
energy between one magnetic domain and the external 
magnetic field: 

W = -c H -H-M (1) 

where ch is a constant, H is the strength of the external 
magnetic field and M is the magnetization of the do- 
main (the positive direction both for M and H is taken 
upwards). Taking into account that the (i — 1) and i 
neighboring domains are oppositely oriented, their total 
energy of interaction with the external magnetic field is: 

W(i) = Wi-x +Wt = -c H ■ H ■ AM (2) 

The quantity AM (the sum of magnetizations of the 
neighboring (i — 1) and i domains) is related to the two 
domains' length's difference (Aa: = Xi — Xi-i) as 

AM = (— 1)S • Ax (3) 

where 7 is constant relating the size of the domain with 
it's magnetization. From (2) and (3) it results: 

W(i) = -chI-H-Ax (4) 

The force Fh acting on domain wall i can be determined 
considering the 5L elementary work performed by this 
force, when the wall is displaced by a distance dl 

SL = F H -dl = -dW(i) = (-lfcHj ■ H ■ d(Ax) = 

= {-ifcHl -H-2dl (5) 

F H = 5 -§ = (-l) l 2c Hl ■ H =(-lT(3-H (6) 

with /? = 2chJ another constant. In our model for the 
sake of simplicity we define the units such that j3 — 1 . For 
positive values of the external magnetic field this force en- 
courages the increase of the domains oriented in the + di- 
rection, and for negative values of the external magnetic 
field this force tends to increase the size of the domains 
oriented in the - direction. 

A second type of force, F m , acting on both sides of 
the domain walls, is due to the magnetic self-energy of 
each domain. This force tends to minimize the length of 
each domain. It can be immediately shown that F. m is 
proportional with the length of the considered domain. 
The Ei magnetic self-energy of a magnetic domain i has 
the form 

Ei = c m M, 2 (7) 
where c m is a constant. As in the previous case 

dEi = c m ■ d{M t 2 ) = -SL = -F m dl (8) 
and from here the F m force: 
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The constant / m is an important coupling parameter in 
this model and acts as the elastic constant of some me- 
chanical springs. 

The system of the F p , F m and Fh forces can be now 
easily mapped on a mechanical spring-block model. 
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FIG. 1. Sketch of the mechanical spring-block model 

The main constituents in this mechanical model arc 
randomly distributed pinning centers, rigid walls sitting 
on pinning centers (describing Bloch-walls) separating + 
and - oriented domains and springs between the walls 
(describing the F m forces). The strength of the pinning 
centers (pinning forces), F p , are randomly distributed fol- 
lowing a normal distribution. Walls can be only on pin- 
ning centers and two walls are not allowed to occupy the 
same pinning center. This constraint implies that the 
number of magnetic domains and domain walls are kept 
constant and are thus a-priori fixed. Domains cannot to- 
tally disappear and new domains cannot appear during 
magnetization phenomena. The clastic springs are ideal 
with zero equilibrium length and with the tension linearly 
proportional with their length. The tension in the elastic 
springs will reproduce the F m forces. Beside the pinning 
forces and the tensions in the springs there is an extra 
force acting on each wall. The strength of this force is 
proportional with the exterior magnetic fields intensity, 
it is the same for all walls but its direction is inverse for 
+ |— and — 1+ walls. This force will reproduce the Fh 
forces. 

The dynamics of this model is aimed to reproduce real 
magnetization phenomena. First N p pinning centers are 
randomly distributed on a fixed length (L) interval, and 
their strengths are assigned. Than a fixed N w number 
of walls are randomly spread over the pinning centers 
(N w <C N p ) and connected by ideal springs. Neighboring 
domains are assigned opposite magnetic orientation. The 
external Fh force is first chosen zero (corresponding to 
H = 0), and we let the system relax to an equilibrium 
configuration. To achieve this we calculate the resultant 
of the F m forces on each wall. If the strength of the pin- 
ning force acting on one wall is smaller than the force 
resulting from the tension of the springs attached to it, 



the wall will jump in the direction of the resultant force 
on the next pinning site. However, if this pinning site 
is already occupied, the wall remains in its original posi- 
tion. We assume that the time needed for the system to 
achieve equilibrium is zero. It is important to note that 
one event (jump) can trigger many other events leading 
to avalanche-like processes. The above dynamics is con- 
tinued until the equilibrium is satisfied for each domain 
wall. The order in which the position of the walls is up- 
dated is random. 

Once the initial equilibrium configuration is reached we 
begin to simulate the magnetization phenomenon. The 
value of the Fh external force is increased step-by-step 
(corresponding to an increasing H magnetic field inten- 
sity) , and for each new Fh value an equilibrium position 
of the system is searched. In equilibrium the magnitude 
of the resulting force on each wall should be less than 
the pinning force acting on that wall. In each equilib- 
rium configuration we calculate the total magnetization 
of the system as: 



(10) 



where U is the length of domain i, and Si is it's orien- 
tation: + 1 for positive orientation, and —1 for negative 
orientation. We increase Fh until no more walls can 
move and the magnetization reaches its maximal value. 
Starting from this we decrease step-by-step the value of 
Fh, and for each new value the equilibrium configuration 
is again reached and the total magnetization computed. 
The system is driven until oppositely oriented saturation. 
From here we increase again the value of Fh and many 
hysteresis cycles are simulated. Throughout the whole 
simulation we assume that equilibrium is instanteneously 
achieved for any value of the H external magnetic field. 

During the simulation we are monitoring the varia- 
tion of the magnetization focusing on the shape of the 
hysteresis loop, jump size distribution, power spectrum, 
Barkhausen signal duration and signal area distribution 
functions. 

The hysteresis loop is the history-dependent relation 
between the magnetization M and the external magnetic 
field H when the value of H is increased and decreased 
successively. The jump size distribution is the distribu- 
tion function for the obtained values of abrupt jumps in 
M throughout many hysteresis loops. The Barkhausen 
signal is given by ^ and it is proportional with an elec- 
tric voltage that would be induced in a detecting electro- 
magnetic coil. Since there is no real time in our simula- 
tions, and the evolution of the system is related solely to 
the driving rate (dH) of the external magnetic field, we 
consider a Barkhausen jump as that is the change 
(jump) in the total magnetization when the driving field 
changes by one dH step. 

We determined the power spectrum of the obtained 
Barkhausen signal by using a Fast Fourier Transform 
(FFT) . Since in our simulations there is no real time, the 
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frequency is also not realistic, and it is defined solely by 
dH which can be considered in a first approximation to 
be constant in time. Thus, we can emphasize in advance 
that the power spectrum determined from this simulation 
should not be considered relevant. 

We also study the shape of the histograms for 
Barkhausen signal duration and signal area. In terms of 
our simulation the signal duration measures the number 
of consecutive dH steps when Barkhausen jumps occur 
(AM/dH is nonzero). Signal area is also related to this 
quantity: it represents the area under AM/dH versus H 
for a nonzero AM/dH sequence. 

The parameters of the model are: N p - the num- 
ber of pinning centers; N w - the number of Bloch-walls 
(N w -C N p ); fm - the coupling constant between the 
neighboring domain walls (corresponds to the elastic con- 
stant in the case of coupled springs); and the dH - driv- 
ing rate of the external magnetic field (change in H for 
one simulation step). The total length of the magnetic 
domains is considered to be of unit length (L = 1) and 
the distribution function for the strength of the pinning 
forces was considered to be normally distributed and this 
distribution has been fixed. We use rigid boundary condi- 
tions: the first Bloch-wall compulsory occupies the first 
pinning center, the last wall occupies the last pinning 
center, and these bounding walls cannot move. This 
constraint means that the geometrical size of our model 
system doesn't change during the simulation. As we al- 
ready mentioned the number of Bloch-walls and mag- 
netic domains are also fixed within this model, domains 
can shrink or grow, but they cannot appear or disappear 
during the simulation. 

III. EFFECTS OF THE FREE PARAMETERS 

Before presenting the simulation results let us discuss 
here the expected effects of the free parameters. We will 
use adimensional units for the forces, their strengths are 
determined relatively to the adimensional strength of the 
pinning forces. The f m , (3, H and x quantities are defined 
through equations (6) and (9) and they are also consid- 
ered to be adimensional. The pinning forces can have 
only positive values and they are normally distributed 
on the [0,1] interval with mean (F p ) = 0.5 and standard 
deviation a = 0.1. 

a. Influence of N p . As N p increases the pinning cen- 
ters are closer to each other which causes many small 
Barkhausen jumps. Concurrently, small number of pin- 
ning centers result less but bigger jumps in magneti- 
zation. From the above arguments one can conclude 
that the value of the parameter N p influences directly 
the shape of the hysteresis loops and the obtained jump 
size distribution histogram. As N p increases the simula- 
tion is more and more time consuming, since many small 
jumps and thus many intermediate equilibrium positions 
are possible. This is thus a first factor that limits the 



value of the used N p parameter. We have performed our 
calculations with N p between 1000 and 5000. 

b. Influence of N w . N w determines the number of 
magnetic domains in our model system. Large values 
of N w require long computation times since the equilib- 
rium becomes more and more sophisticated. Because N w 
should be much smaller than N p , increasing N w would 
lead also to large N p values, which again makes the sim- 
ulation technically difficult. In our simulations we con- 
sidered N w — 100. The parameter N w determines how 
strongly the springs are stressed. Small value of N w 
means that the walls are far from each other and the 
coupling springs between them are strongly tensioncd. 
In this case the obtained avalanches in wall movements 
are usually longer. The N p /N w ratio is one of the most 
relevant quantities regarding the outcome of the simu- 
lations. We consider the N p /N w ratio to be relatively 
"small" if it is between 10-20 and "large" if it is greater 
then 30. For small N p /N w ratio springs are not very 
stressed, thus small number of jumps will occur during 
the magnetization process and strong external magnetic 
field is needed in order to make the walls jump. After 
relatively few jumps saturation is reached and the walls 
will form " pairs" that can be destroyed only by inverting 
the external magnetic field's direction. For large N p /N w 
ratio the existence of the relatively many pinning centers 
causes many small jumps. The small number of Bloch- 
walls causes springs to be stretched and favors the occur- 
rence of jumps even for weak external magnetic fields. A 
large number of steps is needed until saturation is reached 
and the walls are stopped by their neighbors. In simu- 
lations we varied this parameter between 10 and 50 and 
studied it's influence on the BN statistics. 

c. Influence of f m . Since this parameter acts in our 
spring-block model like the elastic constant of the springs, 
its value determines the value of attractive forces between 
neighboring walls. As f m increases the coupling becomes 
stronger, and weak external magnetic field is enough to 
make the walls jump. For small f m values the springs 
are weakly coupled, so a stronger external magnetic field 
is needed to make the walls jump. In our model the 
N p /N w ratio and the f m parameters are strongly related 
to each other. In the case of many pinning centers and 
relatively small number of walls (equivalent with N p /N w 
large) even for weak coupling (f m around 10) many jumps 
occur and equilibrium states are easily reached. When 
the N p /N w ratio is small (around 10) weak coupling re- 
quires strong external magnetic fields in order to make 
walls jump and only a relatively small number of jumps 
are possible. In this parameter region the hysteresis loops 
have only a very limited number of jumps and these 
jumps occur only for high H magnetic field values. If 
the coupling gets stronger (still N p /N w low) equilibrium 
states are very difficult to reach, walls jump back and 
forth and hysteresis loops are totally damaged. When the 
N p /N w ratio is large (around 30), this means that there 
are many pinning centers and relatively small number of 
walls. In this case weak coupling (f m around 10) will be 



4 



enough to make walls jump even for low H values because 
the springs are strongly stressed. The expected result is 
the existence of many small jumps along the whole hys- 
teresis loop. For stronger couplings equilibrium is again 
difficult to reach and the hysteresis curves are expected 
to be damaged. It is obvious that finding the optimal 
parameter region for the simulation is crucial. Wrong 
parameters that generate damaged hysteresis loops, or 
situations where the equilibrium states are difficult to 
reach (walls jump back and forth many times consecu- 
tively) will cause significant artifacts in the jump-size, 
signal-duration or signal-area distribution functions. 

e. Influence of the driving rate in H (dH). This pa- 
rameter is also important for all the considered distribu- 
tion functions. It is obvious that small dH (around 0.001) 
steps produce very many "short" and small jumps, while 
larger steps (dH between 0.01-0.005) make possible only 
the "longer" Barkhausen-signals. The latter means that 
larger dH steps lead usually to many consecutive jumps 
(i.e. large avalanches). Very small dH value will divide 
the bigger and longer jumps into many tinny jumps. In 
the experiments greater dH steps correspond to the case 
where the ferromagnetic sample is submitted to a rela- 
tively fast changing driving field, thus it reaches the sat- 
uration magnetization after few number of large jumps. 
In the experiments it is reported however [2,4,8,15] that 
they used very low frequencies (< 1 Hz) for the driving 
field. This is the reason why we have chosen to make 
the simulations with relatively small dH steps (0.001) 
that corresponds to quasi-stationary driving and allows 
the system to relax to a closer equilibrium configuration 
during each step of the magnetization - demagnetization 
process. 

Pondering the effects of all the above described pa- 
rameters, we have found that the best parameter region 
that produces simulation results in agreement with the 



experimental ones is the following: N p = 1000 — 5000, 
N w = 100 (N p /N w = 10 - 50), f m = 10, dH = 0.001. 
The simulation results presented in this paper are all ob- 
tained for this parameter region. The relevant distribu- 
tion functions were obtained by averaging on 10 indepen- 
dent configurations (different initial states for the walls 
and pinning centers) and considering 10 hysteresis loops 
for every configuration. 



IV. SIMULATION RESULTS 

Characteristic hysteresis loops are plotted on Fig. 2. 
The shape of the obtained hysteresis curves satisfies our 
expectations and fulfills all the requirements for real mag- 
netization phenomena. On these curves one can detect 
many discrete jumps with different sizes, thus the model 
exhibits BN. In addition, when the sample is driven con- 
secutively through many hysteresis cycles the magnetiza- 
tion curves do not follow exactly the same path, although 
the parameters of the simulation were unchanged. The 



qualitative shape of the hysteresis curve is quite stable 
for a wide range of the free parameters. 

The jump size distribution histogram in our simula- 
tions corresponds to the avalanche size distribution his- 
togram in experiments and it is the most relevant distri- 
bution for the characterization of the Barkhausen noise. 
For N w = 100, 7V p =1000, 3000 and 5000, f m = 10 and 
dH = 0.001, the jump-size distribution function 



/(*) = 



N(x, x + dx) 
N t ■ dx 



(11) 



is plotted on Fig. 3. (N(x,x + dx) is the number of mag- 
netization jumps with sizes between x and x + dx, and N t 
is the total number of jumps during the measurement). 




FIG. 2. Hysteresis loops obtained by simulation with pa- 
rameter values: N p = 3000, N w — 100, f m = 10 and 
dH = 0.001. 



• • Np-1000 
b- -a Np=3000 
A- a Np=5000 
— Slope-1.11 




10 100 
x (arbitrary units) 



FIG. 3. Jump size distribution function for N p =1000, 3000 
and 5000, N w = 100, f m = 10, dH = 0.001. The solid line 
indicates a power- function with slope = -1.11. 

Based on this graph we can conclude: (i) for small val- 
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ues of N p the histograms does not show a clear power-law 
distribution for the jump sizes; (ii) the curves for increas- 
ing N p numbers suggest however that for larger system 
sizes the jump size distribution approaches a power-law; 
(iii) taking into account this " straightening" tendency for 
increasing N p values we have fitted a power-law function 
with slope = —1.11 as a guideline. The region where 
this scaling is valid extends however to only one decade 
in our simulations and represents the trend that is to be 
followed if N p is increased to infinity; (iv) the cutoffs at 
the ends of the histograms seem to be finite-size effects. 

The fact that our simulation results indicate a power- 
law distribution in the thermodynamic limit means that 
jump (avalanche) sizes show scale-invariance, a feature 
that is expected to be common for all types of crackling 
noises. 




1 10 100 1000 

f (units in 1/dH) 



FIG. 4. Power spectrum of the simulated BN (N p = 3000; 
N w = 100; f m = 10). 

The power spectrum obtained from simulation suggests 
a white noise over two decades (Fig. 4.). We emphasize 
it again that this is not a real power spectrum since we 
do not have real time in simulations, and thus we cannot 
define frequency. Time evolution is substituted with the 
driving rate dH of the external field and it is considered 
that equilibrium is reached for each simulation step. We 
have plotted thus the "power spectrum" in terms of the 
1 / dH - type " frequency" , assuming the unit time as the 
time needed to change the external magnetic field by dH. 

On Fig. 5. we plotted the signal duration distribution 
functions for different system sizes (1000 — 5000). As 
it results from the graph in our simulations the longest 
jumps lasted around 10 dH steps, we had thus only one 
decade of data. However, the tendency is obvious again: 
the larger the system is, the more the distribution func- 
tion tends to reach a power-law. We have plotted as a 
guideline the power- function with slopc= —2.71 that fits 
the central part of the N p — 5000 histogram. The reason 
that only relatively "short" jumps occurred in our simu- 
lation is that the driving was chosen to be quasi-static. 



The value of dH was chosen small (0.001) and thus the 
jumps arc small and short. Most of the time one dH step 
is not enough to provide the needed amount of energy for 
a further jump if in the previous step a jump occurred, 
thus the big majority of Barkhausen jumps lasted only 
for one dH unit. 




x (arbitrary units) 

FIG. 5. Distribution function for 

signal duration (N w = 100; f m = 10; dH = 0.001). The 
solid line indicates a power-function with slope= —2.71 

On Fig. 6. we plotted the signal area distribution func- 
tion for increasing system sizes (N p = 1000 — 5000). The 
curves suggest a clear tendency: the larger the simulated 
system is, the better a power-law approximation is. From 
our simulation data we would predict for large systems 
a signal area distribution function which is a power-law 
with exponent around —1.15. The cutoffs in our simula- 
tions are due to the finite size effects. 



1.00000 




1000 
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FIG. 6. Signal area distribution function. For the specified 
N p values we used N w = 100; f m = 10; dH = 0.001. The 
solid line is a power-function with slopes —1.15. 
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V. DISCUSSION 

Comparison of the obtained simulation results with the 
experimental ones yields the following conclusions: 

1. The simulated hysteresis loops are visually in good 
agreement with the ones obtained in experiments. The 
best curves are obtained with the parameters: N p = 
5000, N w = 100, (N p /N w = 50), f m = 10 and dH = 
0.001. On the simulated hysteresis loops one can observe 
many Barkhausen jumps with various sizes, just as it is 
expected from the experimental results. 

2. The shape of the simulated jump size distribution 
function predicts in the thermodynamic limit (N w — > oo, 
N p — > oo, N p /N w — > oo) a tendency towards a power- 
law. Experimental avalanche size distribution functions 
exhibit a power-law behavior with an exponent around 

— 1 to —1.3 [14,15]. Our results on relatively small sys- 
tems suggests a power-law behavior with an exponent 

— 1.11. This result is in a reasonable agreement with ex- 
periments and proves that our model is able to account 
for the most important statistics of the BN, although the 
simulated system's sizes were small and finite size effects 
were strong. 

3. As already emphasized, the power spectrum ob- 
tained within our model is not relevant, since we do 
not have real time in our simulations. We considered 
that the equilibration of the system takes place instan- 
taneously and the detection of the Barkhausen signal is 
without any inertia. However, our result which suggests 
a white noise (exponent of the power spectrum w 0) is 
in agreement with O. Narayan's [12] prediction for one- 
dimensional systems. Comparing our results with the 
experiments we found that experiments performed with 
a SQUID device [4] yield also a power spectrum closer 
to white noise than l/f noise, which is thus in agree- 
ment with the present simulation results. The very dif- 
ferent experimental results (power spectra showing white 
noise, l/f noise and even Brownian noise) allow us to 
conclude that the shape of the power spectrum depends 
on the experimental setup and it is not relevant to the 
Barkhausen noise. Similar conclusions were suggested in 
[16]. It seems that the shape of the experimentally mea- 
sured power spectra depends on the inertia of the pickup 
coil. This may be the reason why experiments performed 
with SQUID device [4] indicated that BN is closer to a 
white noise than l/f noise. SQUID devices can follow 
quickly, without any inertia the changes in the magneti- 
zation and this corresponds to the dynamics involved in 
our simulations. 

4. The experimental results for the signal duration 
and the signal area [3,8] distribution indicated a power- 
law behavior for both quantities. For the signal dura- 
tion distribution the measured exponents were —2.2 [3] 
and —1.64 to —1.82 [8]. Although there is a strong dif- 
ference between these results, the experiments agree in 
the validity of the power-law distribution. Our model 
also suggests a power-law distribution in the limit of in- 



finite systems and the predicted exponent is w —2.71. 
For signal area distribution experiments obtained power- 
laws with exponents from —1.7 to —1.8 [3] and —1.74 to 

— 1.88 [8], values that are in good agreement with each 
other. Our simulations predict however a value around 

— 1.15. The numerical agreement is thus not too good, 
but the power-law tendency suggested in our simulations 
can explain at least qualitatively the statistics. 

VI. CONCLUSIONS 

As a conclusion, we can affirm that the model pre- 
sented in this paper is able to describe and explain quali- 
tatively the characteristic features of the BN. The model 
is realistic and reproduces in a pedagogically simple man- 
ner the microscopic dynamics of the magnetic domain 
walls. For a quite broad parameter region our simula- 
tion results proved to be at least in qualitative agreement 
with the known experimental results on the statistics of 
the BN. Despite of the encouraging results the model is 
not perfect. One very important feature is the absence of 
the temperature as parameter. Also, there is no real time 
in simulations, and only the value of the dH magnetiza- 
tion step determines the rate at which time evolves in our 
simulations. The model doesn't account for the pinning 
mechanism and the strength of the pinning forces. It is 
an oversimplified one-dimensional approximation for the 
complex three-dimensional domain topology. Seemingly 
the most serious problem of the present model is that 
the number of domain walls is a priori fixed and domains 
cannot appear or disappear during the dynamics. De- 
spite all these deficiencies this model offers a simple and 
visual picture of magnetization phenomena reproducing 
qualitatively well the statistics of BN. 
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